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MBTtiOD Or CALfBRATJhJG SL4S DRIFT V?fiTK 7EMPES5ATURE FOR 
A VIBRATING STRUCTURE GYROSCOPE. 

This invention relates to a method of talibrating bias drift with operab'ng 
fempersfhjfs over an opening tsfnpaalum range ifor a vibrating stwcture 
. s ' syroscope. 

Vibrating stmcture gyroscopes are known in fnany forms such as tuning 
foric, cylinder or planar ring structures. Planar ring structure gyroscopes can be 
made of silicon with a siiicon vibrating ring structure movably mounted by a 
number of radlaily corr^fiant legs on a fixed support. Drive means are provided 
10 for putting and malnteining the ring structure in vibratory motion at the natural 
resonant mode of the stRicture. Pick-oif mean$ are induded for detecting 
vibration of the vibrating etnjclure, giving a signal proportional to velocity. 

Such a silicon planer ring structune gyroscope is relatively elmpie to 
manufacture in bulic with low production costs but only has a limited bias drift 

15 ~ peribrmance of alsout 1 degree per second over its typical operating 
temperature range of between -40°C and +80°C. Whilst this performance is 
satisfactory for dv|| automotive maricet needs it would be advantageous if the 
bias drift performance could be improved to the order of between 1 to 10 
d^raes per hour whioh would enable such an Improved performance 

20 gyroscope to be utilteed for other applications such as navigation. Bias drift Is 
defined as the variation of the gyroscope output signal as a function of time and 
temperature in tine absence of an appUed rotation rate. 

A simple way of improving the performance of the bias drift over 
temperature would be to model the bias with an external temperature sensoi" 
25 with the use of a polynomial representation as below: 

;&=s«o+tfj5'(5D+tfj5(27*+^?,^(r)»+.„ ---(1) 

where S(T) is the signal representing templerature finom tfie temperature sensor 
and 6 the bias. In order to £^ny out calibration 'process 1h9 'i'lbrating 
30 structure gyroscope is placed Inside a temperature chamber and over a period 
of time cycled over the operational temperature envelope, data-logging both the 
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faias (while the gyro is sEsttonary) and the output of the temperature ssnsor. The 
rate sensor !s held fixed during the period of tesSng. 

A multiple linear regression (?uoh as available from typical spreadsheet 
such as Excel) can 00 used to detsrmlns the cssfficferrts of the m a l^st 
5 squares fit manner fram data over the operational temperature range. It is then 
. assumed that the fit parameters a^, 32 etc. are constant with timei so that this 
one off calibration can then be used for the tested callbrattan gyroscope In 
operation. 

The problem with this temperature sensor approach is that it is not 
• 10 closely related to the gyroscope oharacterlsflcs which accounts for the Bmited 
scope of Improvement that can be achieved. A temperature sensor will not be 
In exacUy the same physical poslBon as the gyroscope so there will inevitably 
be temperature differences espedaBy under dynamic temperature 
environments. This will result in either a lead or lag of the temperature sensor 
15 with respect to the gyro results In apparerrt hysteresis in the modelling with 
increased errors as shown in FTgure 7 of the accompanying drawings, as will be 
later described 

There is thus a need fbr a generally improved method of calibrating bias 
drift with temperature for a vibrating structure gyroscope which at least 

20 minimises these problems. 

According fo one aspect Of the present invention there is provided a 
method of calibrating bias drift with operating temperature over an operating 
temperature range for a vibrating stmcture gyroscope having a substantially 
planar, substantially ring shaped silicon vibrating staictura, primary drive means 

25 for putting and maintaining the vibrating stmcture In earner mode resonance, 
and secondary drive means for nulling response mode motion of the vibraSng 
structure, which secorulary drive means includes means to separate a detected 
response mode motion signal Into a real component induced by applied rotation 
. of the vibratirig structure gyroscope and ..a quadrature component which te an 

30 error term indicative of error mismatch between earner mode resonance 
frequency and response mode resonance frequenty. Including the steps of 
measuring, over an operating tamperatjre range of Hxq vibrating structure 
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gyroscope, primary drive means voltage P which id a measars of chango in 
quality factor Q of the vibrating structure with temperature, vibrating structure 
frequency f which Is a measure of change of temperature of fha vibrating 
stnjctMfs, secondary drive quadrature component values Sq which is a measure 
8 of real component bias enors wBh temperaaji-e. and secondary drive real 
component values Sr which is a measure of change in bias, that Is the zero 
inertial rat© offset, of the vibrating structure gyroscope with temperature 
substituting the values obtained in the relationship 

10 Where aa« are bias calibration ooefBcfents for the vtoiating structure gyroscope 
over the operating tempertrtur^ range, and calculating from the relationship the 
coefficients a,^ to provide a set of bias calibration coefficients for the vibrating 
stmcture gyroscope overthe tested operating temperature range. 

Preferably the coefficients are calculated from the relationship by 

IS carrying out a multiple linear regression on the relationship. 

Conveniently the cbefRdents aa« are calculated by Kalman filtering. 
For a better undeistandlrig of the present invention, and to show how the 
same may be earned Into effect, reference will now be made, by way of 
example, to the accompanying di-awlngs. In which: 

20 Figure 1 Is a graphic illustration of the nomnat mod© vibration pattern for 

a silicon vibrating stmcture gyn>scope, 

Figure 2 is a graphic illustration of the vibration pattern for a silicon 
vibrating stmcture gyroscope when there is mIsaUgnment between the normal 
modes and the primary drive and pick-off means by an angle oc. 

25 Figure 3 is a graphic illustration of the Corriolls Effect on the primary 

mode vibration pattern for a silicon vibrating structure gyroscope, 
* Figure 4 is a dfagrammatic representation of a "conventional electronic 
control circuit for a conventional vibrating structure gyroscope having a silicon 
ring shaped vibrating structure. 
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Flgure 5 Is a block function diasram of the cDntrol dfcult of Figure 4, 
Figure 6 Is a block ftjncfion diagram of a modrfied version OFf the control 
drcuitof Rgure4. 

Rgur© 7 is a graphical representation of bias versus frequeticy sliowing 
the hysteresis effect on bias over temperature. 

Referring now to Figures 1 to 4 of the accompanying drawings a silicon 
ring shaped vibrating slmclure 1 has two lUndamental modes, differing by. 
typically, less than 1 Hz in frequency. Ideally the primary mods is driven at tts 
natural ftequency, and the secondary mode is suppressed- In practice the ring 
shaped structure 1 oscillatBs at a frequency between the primary end 
secondary natural fiequencies, and both modes are excited. The magnitude of 
oscillation is held approximafely constant across the full operating temperature 
range. An AGC(automatlc gain control) loop 2 operates to control the amplitude 
of the primary moiton. 

The gyro operating principle involves the use of Coriolis tbrces coupling 
the primary and secondary modes when inertial rotation occurs about an axis 
normal to the plane of the ring shaped stmcture 1. The force needed to 
suppress the Coriolis Ibnse is proportional to the inertial angular rate. 

The main enor driver Is the angular misalignment of the primary drive 
and the primary mode as shown in Figure 1. This misal^nment angle varies 
acn>ss the temperature range, and is subject to hysteresis. It Is not simply a 
function of present temperature, but is a hysteretic, stochastic tjnction of 
temperature history, which Is unknown, it Is this hysteresis that gives rise to 
errors that cannot be calibrated out off-line by conventional means. Apart from 
the output of measured gyro rate 3 (conupted by unknown errors), additional 
measurements are available from the primary drive voltage 4, frequency 5, and 
secondary quadr^ure drive or bias 6. These measurements are known to be a 
function of the various error drivers. 

Figure 1 of the accompanying drawings shows the primary normal mode, 
30 with primary drive means 7 and primary pick-off means 8. Also shown, at 45 
degrees to the primary mode, are the axes of the secondary mode with 



^miiji^i 



lCW-2002 16:33 • FgM Ir,t^l lectual R-op pf^^ oFFtCE 



P. 10 



secondary drive means S and secoRdary pick-off means 10. The mode posltbn 
and driver and pick-off posltioris are all shown fn their ideal positions. The 
secondary mode Is rotated 46 degrees relative to the primary, 

MisaHgn?nsnt of t«iB normal nmodes wffti the drivers and pick-offs by an 
5 anglea, as shown In Bgurs 2. results in . a combination of the primary and 
secondary modes being used to null moflon at 45 degrees to the primary drive 
means 7. The primary pick-off means 8 will delect mainly the prlmaiy mode, 
but also the secondaiy mode, to a lesser degree.. The primary drive is the 
reference direction from the gyro centre, from which all oth^ angles are 
10 measured. 

Due to nominal symmetry of the ring shaped striiciure 1, it appears that 
the prfrnary mode could take up any position on the structure 1, however minor 
asymmetries may determine Its actual position^ Thus the mode position tends to 
be locked to one angle a due to anisotropies around the ring shaped structuns 
15 1. U, In Figure 2 represents the primary mode amplitude and 11 the primary 
mode node. 

* Rgure 3 shows the velocity vectors 2 at four points for the primary mode. 
Due to Coriolls effect in the presence of rotaSon, there are D'Alemberf s forces 
13 applied at each point Clearly, Vne primary mode motion in the presence of 
20 rotation gives rise to forces, which drive the secondary mode, and by symmetry, 
the serandaiy mode moMon in the presence of rotation gives rise to forces, 
which drive the primary mode. 

The pattern of Coriofis acceleration suggests that the secondary mode 
wUI be driven by D'AIembert's forces. This leads to a cross-coupling between . 
25 the primary and secondary normal modes, in both directions. There is a gafn 
associated with this effect, the Bryan gain. Thus, for a rotating ring staicture 
there Is a special effect, referred to here as the Bryan Effect, which ibrms the 
physical ba^is of the gyro's operation. 

Typically a sllicon"fing shaped vibrating structure gyroscope is controlled * . 
M by an electronic control system as shown In Figure 4 of the accompanying 
drawings. This has a primary loop 14 and e secondary loop 15. IVlaking up the 
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primary ioop 14 is a frequency corrtroi Socsp 16 and an amplitude control loop 17. 
The loop 16 contains a phase locked loop 18 and a voltage controlled oscillator 
1 9. The loop 17 contains an automatic gain control 20 to stabiFise the ampHtuds 
of mcScn. hsvlng an Input V/^gc^ voltage to set the amplitude of motion. Loop 14 
5 Is connected between the primary pickKjff means 8 and the primajy drive 
means 7. The secondary loop 15 is connected between the secondary pick-off 
means 10 and the secondary drive means 9 and acts as a 
demodulating/remodulatlng nuD control loop. To ttifs end K contains a 
quadrature component loop 20 and a real component loop 21 . 

10 Figure 5 shows a block function diagram of the control system of Figure 

4. In which like parts have been given like reference numbers. The output from 
the automatfc gain control 20 is indicated at 23. 

Figure 6 shows a block function diagram of the control system of Rgure 
4, In which like parts have been given like reference numbers, but also in which 

15 there la an extra servo loop 24 connecting the phase tocked loop 18 to the 
gyroscope vibrating strut^jre 1 to ensure that there is a SO* phase shift between 
the primary drive voltage 4 and the primary plofc-off signal- 
Over temperature, due to the change of modulus of ttie sili<X3n and the 
size of the ring structure 1. there is a char^ge of resonance frequency of 

20 -0-4HZ/C, With a typical phase error of the loop of 0.5 degrees, the temperature 
of tlie ring structure can be determined to an accuracy of 0.01 Hz corresponding 
to 0.026 C The frequency of vibration of the structure 1 can be determined 
very accurately during calibration and this Implies that the temperature 
environment of the ring structure can be very precisely detemnined with no time 

25 lag between the "temperature sensor^ (the frequency of the ring) and the ring 
structure itself. 

The frequency of the ring structure 1 when It is vibrating at resonance Is 
a good measure of the actual temperature of the ring structure with a typical 
coefficient .of -0.4Hz/C. The frequency, and the gyro bias with th6 gyroscope 
30 stationary In a temperature chamber can be recorded over the operational 
range, which Is typically -40C and an Improved model ibrthe gyro bias can be 
written as: 
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w^ere f Is tfie frequency of the rmg. This again is a csdibration whrch {s carried 
out overtiie operationai temperature range of the gyro wi&i data b&ing coii^eul 
for the gyro bias b, scale factor s and frsquenQr f of the ring structure as. a 
function of time when a typical femperature profDe as In Figure 7 Is found. 

Whilst this model results In an improvement In the gyro bias performance 
over the extemal temperature recording process, it cannot eliminate the effects 
of hystetiesis of bias over ferapenature as thens is no time history In the model. 

Figure 7 shows that for variation in temperature tfiere is a hysteresis 
efTect in that the solid line In the gnaph shows a deviation depending upon 
whether the temperature is increas&ig or decreasing, as shown by the arrows, 
fiT>m a mean track 25. Hence this process fbr calibrating bias drift is inherently 
inaccurate. 

With reference to the block diagram In Figure 6, tiiere are a number of 
measurands in the gyro as follows: 

a) primary drive voltage 4, given *ie symbol P 

b) Frequency 5. given the symbol f 

c) Secondary quadrature drive levels 6, given the symbol Sq, and 

d) Secondary real drive 3, which is the nonrtal rate output Sr 

The primary drive voltage P is a measure of the drive nequlred to maintain Qie 
primary oscitlatfon of vibrating structure 1 at a fixed amplitude. As the 
temperature changes the Q of the gyroscope decreases with increasing 
temperature rsquinng a higher drive level. Thus the primary drive voltage P is a 
measure of 1/Q as a tiinction of temperaturB, and this knowledge can be used 
to take out the effects of temperature variation of Q. The fr^quenty, as 
described before, is a measure of ring strut^ne temperature and represents a 
very accurate indicator. The seoondaiy quadrature signal St, is related to enor 
affecte of the real bias. In particular the ring shaped vibrating structure 1 is 
^own to have two elgen modes . , 



20-NDU-2002 16=34 FROM Intellectual Prop 



TO EB PfiT DFFECE 




5 • 



10 



15 



These modes are not naturally aligned to the primary and secondary 
drive axes, but have a misalignment angle a. When the primary drive is 
energised, if a Is not equal to zero then a linear combination of the two eigen 
modes Is exerted. These aigen modes are characterised by resonant 
frequencies fi and fa and quality factors Qi and Qa. If the two modes do not 
have the same resonant frequency, and the angle a is not equal to zero then a 
quadrature drive signal Is required to null out the quadrature motion. It can be 
shown that the quadrature level Is related to the mode angle a. 

There are also errors in the bias , the zero inerfal rate offeet of the gyro, 
whioh Is measured by Sr- It is known that when the Q factors of the two modes 
are different, then there is a bias error proportional to sin 4a and the values of 
AQ. It Is seen that there is thus a correlation between the value of S<, and the 
proporSfon of fh© error in Sr due to AQ. Thus the variation of S, will track the 
variation of S,, so that the error due to delta Q in gyro bias can be 
rampensated. 

Over the temperature range of ±80' this may be written as: 



where an ts a fit coeffitient liom a least squares fit 

In a similar way the gyro bias will have error temns which vary as the Q of the 
mode changes, which is measured by the primary drive level, and this can be 



' — (4) 



written: 



---(6) 
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Puiiing aS ihee© terms tcg^er we can writs: 

Note that In this expression the values of K I and m can be both positive 
and negative, 

in order to cover the major entMS that are known on the gyn:> typicaily a 
total of Ttenns are used . 
These cover enors such as: 

a) phase enor of the primary drive with resped to the secondary pidc off. 
These are electronic phase errors with a typical value of 0.5 degrees due 
to limited pre amplifier bandwidth and time delays through integrated 
cifcuite. 

b) Effects of delta Q as discussed above 

c) Effects of the angle offeet between the secondary pick off and primary 
drive 

15 d) Effects of bias voltage to null the original bias offset ^ 

e) Eflects of Inductive pick up between primary drive and secondary pick off 

f) Effects of capadtive pick up between primary drive and secondary pick 
off 

g) Effects of magnetic field anomalies. 
By a sullable Choice of the powers of k, I and m It is possible to model 

these 7 error terms using the three measurands. During the calibration process 
the values of the measuiands - Sr Sq P and f are data-logged as the 
temperature varies over a range of ± SO'C. A multiple linear regression is 
carried outio determina tiie oosfilclents aum* 

Other analysis technkjues can be used to derive the coefficients irom the 
data such as Kalman filtering techniques. 



25 
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As fh© measurands are correlated to the gyro bias errors drivers, 
measuremente of the these terms will improve the bias compensation. The 
polynomial ftwm given above is then linked back to physical error sources (as 
discussed), and the closer the calibration is to the physical reality then the 
5 better thamodeiiing that can be achieved. 

As a result it is now Is posstoie to model gyro bias hysiereslSt which a 

t 

simple polynomial model cannot achieve. This Is du© to the correlation of real 
and quadrature bias in particular. Therefore If there is hysteresis in real bias 
there will generally be a corresponding hysteresis in quad bias as the earor 
10 mechanisms (such as variation of alpha and hysteresis of alpha with 
temperature) are common to these two. 

In practica the callbretlon coefficients atim are established for a specific 
gyroscope by applying in an oven, a desired temperature sequence to the 
gyroscope. Test date, incdudlng the pre\flously specified measurands is 
1S recorcled and a set of bias calibration coefRciente established for the gyreiscope. 
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CLAIMS 

1. A me^od of calibrating bias dilft with operating temperature over an 
operating temperature range fbr a vibfHisf^ siTuciure gyrcsospe ha>4ng a 
5 substantiaHy pl^ar. substanfi^y rkig shaped siiicon vibrating stmcture, primary 
drive means for putting and maintaining tits vibrating i^ctUFe in carrier, mode 
resonance, and secondary drive means for nulling response mode motion of tiie 
vibrating structure, whicii secondary drive means indudes means to separate a 
detected response mode motion signal iriio a real component Induced by 
10 applied rotation of ttie vibrating structure gyroscope and a quadrature 
component whicii is an error term indicative of error mismatch between carrier 
mode resonance frequency and response mode resonance frequency, indudtng 
the steps of measuring, over an operating temperature range of the vibrating 
structure gyroscope, primary drive means voltage P which is a measure of 
15 change in quality fector LQ of the vibrating structure v/t&i temperature, vibrating 
structure frequency f which is a measure of change of temperature of the 
vibrating structure, secondary drive quadrature component values which Is a 
measure of real component bias errors with temperature, and secondary drive 
real component values Sr which is a measure of change in bias, that is the zaro 
20 inertia! rate oi^set, of the vibrating stmcture gyroscope with tem(:>erature, 
substituting the values obtaft^ed in the reiationship 



whare ajom are bias calibration coeffldents for the vibrating structure gyroscope 
over the operating temperature range, and calculattng from the relationship the 
25 coefficients aiam to provide a set of bias calibration coefficients for the vibrating 
structure Qyrosoope over the tested operating temperature range, 

2. A method according to claim 1., in which the coefftcients are 
^ calculated from the relationship by carrying out a muftiple linear. regression on;,, 
the relationship. 
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3. A method accorxjlng to claim 1, in which the coefficiente atew are 
calculated by Kalman fJlterifig. 

4. A method of calibrating bias drift with operating temperature over an 
operating temperature range for a vibrating structure gyroscope, substantially 
as hereinbefore described. 



• 
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ABSTRACT 

METHOD OF CALIBRATINe BIAS DRIFT WITH TEMPERATURE FOR A 

VIBRATING STRUCTURE GYROSCOPE 

5 A method of calibrating bias drift with operating temperatuns over an operating 
temperatuVe range for a vibrating structure gyroscope havirsg a substantially 
planar, substantially ring shaped silicon vibrating stmcturs (1), primary drive 
means (7) for putting and maintaining the vibrating structure (1 ) in carrier modo 
resonance, and secondary drive means (9) for nulling response mode motion of 

10 the vibrating strwdure is described. The secondary drive means (9> includes 
means to separate a detected response mode motion signal into a real 
componerrt Induced fay appOed rotation of the vibrating structure gyroscxjpe and 
a quadrature component whiot* is an error temn indicative of en-or mismatch 
between carrier mode resonance tirequenty and response mode resonance 

15 frequency. The method includes the steps of measuring, over an operating 
temperature range of the vibrating structure gyroscope, primary drive means 
, voltage P which te a measure of change In qualHy fector LQ of the vibrating 
structure w/ith tsmperatut«. vibrgtlng stnicture frequency f which is a measure of 
change of temperaftire of the vibrating stmcture. secondary drive quadrature 

20 component values Sq which is a measure of real component bias errors wKh 
temperature, and secondary drive real component values Sr which is a measure 
of change in bias, that is the zero inertlal rate oifeet. of the vIbraUng stnicture 
gyroscope with temperature, substituting the values obtained in the relaUonshIp 

tin 

25 where ate bias calibration coefficients tor the vibrating stmcture gyroscope 
over the operating temperature range, and calculating from the relationship the 
coefficients aon to provide a set Of Was calibration coefficients for the vibrating 
stnjcture gyroscope over the tested operating temperature range. 
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